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Abstract 

The purpose of this paper is to estabhsh that for any compact, connected C°° Riemannian 
manifold there exists a robust family of kernels of increasing smoothness that are well suited for 
interpolation. They generate Lagrange functions that are uniformly bounded and decay away 
from their center at an exponential rate. An immediate corollary is that the corresponding 
Lebesgue constant will be uniformly bounded with a constant whose only dependence on the 
set of data sites is reflected in the mesh ratio, which measures the uniformity of the data. 

The analysis needed for these results was inspired by some fundamental work of Matveev 
where the Sobolev decay of Lagrange functions associated with certain kernels on $7 C R'^ was 
obtained. With a bit more work, one establishes the following: Lebesgue constants associated 
with surface splines and Sobolev splines are uniformly bounded on M.'^ provided the data sites 
S are quasi-uniformly distributed. The non-Euclidean case is more involved as the geometry of 
the underlying surface comes into play. In addition to establishing bounded Lebesgue constants 
in this setting, a "zeros lemma" for compact Riemannian manifolds is established. 



1 Introduction 



Radial Basis Functions (RBFs) and their counterparts on the sphere (SBFs) are now well-known 
and widely used for fitting a surface to scattered data arising from sampling an unknown function 
defined on M'^ ( or S"^). The fitting of the surface is typically implemented by means of interpo- 
lation or discrete least squares. Both of these tools fit under the larger umbrella of kernel based 
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approximation, a methodology, suitable for treating functions defined on very general domains. 
This approach has achieved success in approximating scattered data on spheres and Euclidean 
domains (including the case when data is arranged on lower dimensional manifolds), as well as 
more esoteric domains like graphs and Lie groups [30^ [T4l [TBI [20] . The hallmark of this methodol- 
ogy is to use the kernel k to create an approximant, s, by taking linear combinations of the form 
s = ^^=i-A-jK{-^Xj). Evidently the first challenge is to find a "scheme" assigning the coefficients 
and sometimes the "centers" Xj, j = 1...N. 

One such a scheme, interpolation, is known to work quite well and the interpolants even exhibit 
"local" behavior just as in the case of univariate spline interpolation. The reasons for the local 
behavior of spline interpolation, i.e., where changing the data locally only significantly alters the in- 
terpolant locally, are well understood. In fact since univariate splines are compactly supported, the 
corresponding interpolation (or collocation) matrices are banded. Thus, whenever the interpolation 
matrices are invertible, the inverse matrices decay exponentially away from the main diagonal (see 
[6] and [?])• This fact, together with the local support of spline functions, implies the exponentially 
fast decay of Lagrange functions away from their "centers" . 

In contrast, many popular kernels (including RBFs and SBFs) are typically globally supported 
and the interpolation matrices, 6= := (k(^, C))(g,(;)eH! full. Even in the case of the compactly 
supported Wendland functions, as the number of data sites increase while still interpolating with 
shifts of a fixed 0, the bands of the interpolation matrices become ever larger. 

Nevertheless, in this paper, we will show that for any compact, complete Riemannian manifold there 
exists a robust family of kernels of increasing smoothness, {Km,M : m € Z, m > d/2} (see Section 13.31 
for a definition), generating Lagrange (or fundamental, or cardinal) functions - the functions = 
X^ceH ^C'^('' C) satisfying X£.{C) = ^{C^O ^ that are uniformly bounded and decay away from their 
center at an exponential type rate. An immediate corollary is that the corresponding Lebesgue 
constant will be uniformly bounded with a constant whose only dependence on H is reflected in 
the mesh ratio, which measures the uniformity of the data, a simple parameter measuring the 
"badness" of the geometry of the data (measured as the ratio of mesh norm to separation radius - 
see definitions below). 

The Lebesgue constant, L := L{E) := sup^^^'^^^^\x^{x)\, is the Loo Loo operator norm 
of the projector that maps continuous functions to interpolants, / i-)- /=/ := fiOx^- This 

constant provides a measure of the stability of the interpolation process. This, it should be noted, is 
independent of the myriad ways the interpolant may be determined; how coefficients are determined 
depends strongly on the basis for the space of interpolants - using the kernels as the 

basis for this space may lead to coefficients of indeterminant size, (c^)gGH •= (/(^))^g„, while 
using the basis of Lagrange functions always involves a stable selection of coefficients. Thus, the 
boundedness of the Lebesgue constant indicates the stability of the kernel based interpolation we 
consider, while the rapid decay of the Lagrange functions shows that this kind of interpolation is 
local in nature. 

The challenge of finding stable interpolation processes is underlined by the fact that, for many 
of the basic tools in approximation theory (e.g., polynomials, trigonometric functions, spherical 
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harmonics), interpolation is known to be unstable. A consequence of our main result is that for 
any S'^, the family of SBFs gd give rise to bounded Lebesgue constants (each corresponding, 
roughly, to the fundamental solution of an elliptic differential operator). This result stands in stark 
contrast with any polynomial interpolation scheme on the sphere. In particular, let IIx, be the space 
of all spherical harmonics of total degree L or less and let T be any bounded linear projector (i.e., 
satisfying = T) 

T : C{S'^) Ul 

with both spaces endowed with the L°°(S'^) norm. Then, as shown in [26], 

||T|Ioo > L"^, d>2. 

Thus, if H has mesh norm h and h ^ q, then #H = (D{h~"') w dimll/j-i. In this case any polynomial 
interpolation scheme would have Lebesgue constant > 0((^)'''^~^^^^) • 

Another happy consequence of the uniform boundedness of the Lebesgue constant is the so-called 
Lebesgue lemma, which states that interpolation is near best L^o approximation from the space of 
interpolants, 5= = span{K(-,^) : ^ G H}. Since /=s = s for each s € 5=, one has 

11/ - /h/II <\\f-s- Hf - s)\\oo < (1 + i)dist(/, 5h). 

This is significant for kernel based approximation, because L^o error estimates for interpolation 
have been notably imprecise. By using direct methods of error analysis, it has, to date, only been 
possible to measure error for target functions coming from a certain reproducing kernel Hilbert 
space associated with k. (Frequently this space is an L2 Sobolev space; for an example of this kind 
of error estimate, see Lemma [3.91 of this paper.) Because one cannot assume that data is generated 
by a target function from a specific smoothness space, it is desirable to have an accurate view of 
approximation power for target functions of arbitrary smoothness, e.g., by having error estimates 
for a scale of spaces capturing fractional smoothness. Unfortunately, precise error estimates for 
interpolation from the scale of Holder spaces (C*), Bessel potential spaces (W^), Besov spaces 
(-^00, g) other spaces measuring smoothness in Lqo have not been forthcoming. Recently, this has 
been addressed for kernel based approximation, using schemes other than interpolation, [17\ I24j. 
and the results of this article point the way toward fixing a significant gap in the study of kernel 
interpolation. 

To illustrate this, we compute the SBF, K2§2, on the sphere in M'^: S^. This is a kernel having 
a bounded Lebesgue constant, but it also belongs to the scheme of kernels considered in [21], 
and, hence, its best approximation properties are well known. It is associated with a fourth order 
elliptic partial differential operator on S^, and provides approximation order a for functions having 
smoothness (J in Lqo'; provided cr ^ 4. This means that for a function, y, having smoothness o", 
dist(/, 5"=) = 0{h"), where h is the mesh norm of H. Previously, the most one could have said about 
approximation power of interpolation with this kernel is that — /||oo = 0{h^) for functions, /, 
having two derivatives in L2. Thus, the Lebesgue lemma ensures that interpolants provide optimal 
Loo approximation orders, and not only for continuous target functions of low smoothness; it also 
significantly improves the approximation power for functions having greater smoothness. 
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1.1 Methodology 



There is a unifying thread that explains the fast decay of the Lagrange functions for both the 
univariate splines and kernels we consider. Both interpolation problems can be considered as 
minimization problems with interpolatory constraints: /=/ = argmin{||s|| : = This 
variational aspect of the theory plays a well-known role in the etymology of the term spline - the 
draftsman's spline is a flexible implement that passes through points fixed on drawingboard by 
assuming a smooth curve that minimizes a "bending energy," the L2 norm of the second derivative. 
The connection between kernel interpolation and variational problems constrained by interpolatory 
conditions in Hilbert spaces is the basis of the approach taken by Madych and Nelson (cf. |21j). 
which has since become a fundamental part of RBF and SBF theory. However, the idea of cleverly 
exploiting this variational property to obtain decay of Lagrange functions was developed by Matveev 
in an impressive paper [23j, that has gone mostly ignored (Johnson's use of Matveev's results to 
achieve surface spline Lp saturation orders, [19], is a notable exception). 

The subject of Matveev's interpolation problem is the D'^-spline, the interpolant minimizing a 

/ \ 1/2 

Sobolev seminorm |/|m,Q = (Sjci|=m/n f{x)'\^dx\ . For general 0, the minimizing kernel is 

computationally intractable, but the case = M'^ gives one of the pre-eminent families of RBFs, 
the well-known surface (or polyharmonic) splines considered by Duchon jlH [TO] and many others. 

Matveev's technique establishes a 'Sobolev decay' of the Lagrange functions, i.e., the decay of 
the Sobolev seminorm of Xi — the Lagrange function for surface spline interpolation — over the 
punctured plane, W^\B{^, R), having a hole of radius R (cf. Section d]). Johnson applied a Sobolev 
embedding argument to obtain pointwise estimates [191 Corollary 2.3] for Lagrange functions whose 
centers come from a fairly regular set of data sites H. It is clear that exponential decay such 
as this holding for all Lagrange functions would be sufficient to uniformly bound the 'boundary- 
free' Lebesgue constant. A minor, albeit oblique, modification of his argument would show this. 
We remark that Johnson's paper had a completely different focus than stability of surface spline 
interpolation, and exponential bounds on the Lagrange functions were ancillary to its main goal. 

In this article we have improved the approach to estimating Lagrange functions. Although it has 
been adapted to handle interpolation on compact Riemannian manifolds, the setting for the core of 
the argument treating Lagrange functions is the tangent space at ^; it is essentially still an argument 
set in R"^. The article is roughly divided into a component where the geometry of the manifold 
is important, where a powerful relationship between Sobolev norms on the manifold and on is 
established, and a second component, in which the Riemannian geometry is not so important, that 
gives Sobolev and pointwise decay of the Lagrange functions modeled on Matveev's argument — 
although with a different underlying variational problem, and with a different embedding argument, 
providing stronger pointwise estimates. 
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1.2 Outline 



Section [2] is devoted to Riemannian manifolds. Since this paper is about approximation theory, 
not differential geometry per se, we provide, in Section 12.11 the requisite geometric background on 
Riemannian manifolds - discussing geodesies, Christoffel symbols, exponential maps, and so forth. 
Section [2?2] concerns the covariant derivative, which is essential for defining invariant Sobolev spaces. 
In particular, we provide various pointwise bounds on such derivatives in terms of ordinary (flat 
space) derivatives. These bounds are what allow us to use Euclidean-space Sobolev estimates to 
obtain similar bounds for the covariant Sobolev spaces on a Riemannian manifold. 

The covariant Sobolev spaces we will work with are due to Aubey (cf. [2]); we define and discuss 
them in Section [3l In Section 13.11 we use the covariant-derivative bounds derived in Section 12.21 
to provide a very strong equivalence between Sobolev norms over regions on the manifold and 
corresponding regions in the tangent space, with constants of equivalence that are independent of 
the geometry of the region. Section 13.21 treats the problem of estimating functions having many 
zeros in a star-shaped domain in domain in W^, and from that and the strong equivalence above, 
deduce similar results for a geodesic ball on a manifold. This is a key element to estimating Lagrange 
functions, which, of course, have many zeros. A consequence of this "zeros lemma" is a conventional 
uniform estimate on the interpolation error for many kinds of kernels. In Section [3. 3^ we introduce 
the family of kernels associated with bounded Lebesgue constants. These are reproducing kernels 
for the invariant Sobolev spaces Ty|"(M) mentioned above. 

In Section we provide pointwise estimates of Lagrange functions and the main theorem, which 
estimates the size of the Lebesgue constant for interpolation on manifolds or on M*^ with surface 
splines or Sobolev splines. In Section [4. 1^ a bound on the size (in Loo) of Lagrange functions from 
a very general class of kernels is presented (kernels having a Sobolev space as their native space) . 
Section I4.2l treats the decay of the Lagrange functions from the class of kernels associated with the 
inner product of Section [3l This is done in roughly three stages. The first stage provides the basis 
for a "bulk chasing" argument by showing that the bulk of the tail of the Sobolev norm is contained 
in a thin annulus near the beginning of the tail. The second stage iterates this result, showing that 
the tail of the Sobolev norm decays exponentially. The third stage uses the zeros lemma to provide 
the exponentially decaying pointwise estimates of the Lagrange function. Section [4. 31 gives the main 
result, which follows in a direct way from the decay of the Lagrange functions. 

The concrete, spherical, example of the SBF K2^s.2 is explicitly calculated in Section [5l It is shown 
that this SBF is a perturbation of the fundamental solution of the bi-Laplacian, on S^, and 
hence is of the type treated in f24j. In particular, it inverts a fourth order differential operator, and 
is, for this reason, not very different from either the thin plate splines in or the univariate cubic 
splines. Since much is known about the abstract approximation properties of these types of SBFs 
(but little is known about the convergence of interpolation), an application of Lebesgue's lemma 
reveals that interpolation provides L^o approximation orders commensurate to the smoothness of 
the function. In particular, 11/=/ — /||oo = 0{h'^) for functions in C^[§'^). 
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1.3 Notation 



M denotes a compact, complete Riemannian manifold of dimension d. More will be said about 
it in Section [2l but for now we simply remark that it is a compact manifold - a Hausdorff space 
possessing a (maximal) family of charts {4'j^Uj)j^g, where each map 4>j : [/j C M'^ — )• V^- C M is a 
bijection and the compatibility condition c/)^"^ o (pj g C°° (^<f>J^ (Uk)) , for every k,j S 3 is satisfied. 
At every point p € M, we denote the tangent space by TpM. This fact gives rise to a notion of arc 
length of curves, a metric: (.^, a) i— )■ dist(^, a), and measure, which we denote by //. 

The basic tools used for estimating Lagrange functions will be the various Sobolev norms and 
seminorms on balls, complements of balls and annuli. 

Although the formulation for Sobolev spaces on manifolds is rather technical, hence postponed until 
the next section, the Euclidean versions are standard and quite simple. For an open set 17 C M'^, 
and a positive integer j, we define the seminorm of the function u : $7 — )• M as 

1/2 



dx 



where the multinomial coefficient is (^) := . The Sobolev norm is then 



1/2 



Frequently we must simultaneously refer to sets on the manifold and on Euclidean space. To avoid 
confusion, we adopt the following notation. We denote the ball in M"^ centered at x having radius 
r by B{x,r). For t > 0, we define the annulus in R*^ centered at x of thickness r and outer radius 
tr by 

A{x, t, r) := B{x, tr)\B{x, {t - l)r). 

We denote the ball in M centered at a having radius r by b(a, r). For t > 0, we define the annulus 
in M centered at a by 

a(a, t, r) := b(a, tr)\b(a, {t — l)r). 

The set S C M will be assumed finite. Its mesh norm (or fill distance) h := sup^.^^ dist(x, H) 
measures the density of S in M, while the separation radius q '■= \ inf^ dist(^, Q) determines the 

spacing of H. The mesh ratio p := h/q measures the uniformity of the distribution of H in M. 



2 Riemannian Manifolds 



Since this paper is primarily about approximation, we will provide the necessary background in 
Riemannian geometry - e.g., Christoffel symbols, covariant and contravariant tensors, exponential 
maps, etc. - that will be needed here. For more information, see the books by do Carmo [8l[9]. 
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2.1 Background 



Throughout our discussion, we will assume that (M, g) is a d-dimensional, connected, C°° Rie- 
mannian manifold without boundary; the Riemannian metric for M is g, which defines an inner 
product gp{-, •) = (•, •)g,p on each tangent space TpM. As usual, a chart is a pair (il, (p) such that 
H C M is open and the map (?!) : H — )■ M'' is a one-to-one homeomorphism. An atlas is a collection 
of charts {(11^, indexed by a such that that M = Uo-ila and, when (pai^a) H 0/3(11/3) / 0, 
0/3 o (p^^ is C°°. In a fixed chart (it, 0), the points p G il are parametrized by p = (p~^{x), where 
X = {x^,. . . ,x'^) G [/ = 0(il). As is common in differential geometry, we will use superscripts to 
denote coordinates. Thus, is the j^^ local coordinate for p relative to the chart. Most of the 
analysis we will do takes place in local coordinates. 

In such coordinates for H, the tangent space TpM at p has a basis comprising the partial derivatives 
(ai?)^' '^'^^^^ i = 1) • • • ) c?- This allows us to define smoothly varying vector fields on U via 

d 

ej(a;) := j = 1, . . . , xeU = 0(il). 

For / G C°°(il), so that / o G C°°{U), the action of a vector field v thus is given by the 
following: 

for any x = (t){p) G U. The coefficients = v^{x) are smooth functions; they are called the 
contravariant components of v relative to the ej basis. We write v(/)(p), rather than v(/)(a;) 
because the expression on the left is independent of the choice of coordinates for p. Obviously, in 
terms of tangent vectors alone we have that 

d 

v = ^t;Jej-. 



Metric tensor Let and be the contravariant components for the vectors v and w in TpM. 
The inner product for the Riemannian metric g is given is then given by 

d 

{v,'w)g,p = ^2 9ij{^)'^^{^)'^'' ^ = ^{p) (contravariant form). 

The gij are the covariant components of g, and are given by gij{x) = {ei,ej)g^p, which are the 
entries of the dx d Gram matrix for the basis {ej}. As was the case with v(/)(p), (v, y^)g,p is itself 
independent of coordinates. In addition, if v and w are C°° vectors fields in p, then (v, w)g,p is 
also C°°. Finally, because gij is a Gram matrix, it is symmetric and positive definite. 

The space dual to TpM is the cotangent space, and it is denoted by T*M. Corresponding to the 
basis Bj = for tangent vectors, we have the dual vectors eP = dxK The two bases are related 
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through the metric via 



&'{x)= ^ g^''{x)ek{x) and ek{x) = ^ gkj{x)&'{x). 
k€{l...d} jG{l-d} 

The matrix with entries g^^ is simply the inverse of the metric tensor (gjk)- Making the usTial 
identification between a vector space equipped with an inner product and its dual space, wc may 
represent a vector v in terms of the e^'s as v = Ylje{i d}'"j^''- ""/^ called 

the covariant components of v. The natural inner product of vectors v = d} '"j^ 

w = J2je{i d} '^3^ expressed in covariant components is 

d 

{v,w)g^p = ^ g^^ {x)vi{x)wj{x) (covariant form). 

An order k covariant tensor T is a real-valued, multilinear function of the fc-fold tensor product of 
TpM. We denote by T^M the covariant tensors of of order k at p. In terms of the local coordinates, 
there is a smoothly varying basis e*i (8) • • • 6**= for the /c-fold tensor product of tangent spaces. 
Thus, the covariant tensor field T of order on it can be written as 

T= Tje'i ® •••®e^^ 

ie{i...d}>' 

where we adopt the convention i = {ii,... ,ik). The Tj are the covariant components of T. The 
metric g is itself an order 2 covariant tensor field. One can also define contravariant tensors and 
tensors of mixed type. 

Because T^M = TpM®- ■ -^TpM (k times), the metric g induces a natural, useful, invariant inner 
product on TpM; in terms of covariant components, it is given by 

(S,T),,p= E 9''''---9''''S,T^. (2.1) 
i,ie{i...rf}'= 

The corresponding norm will be denoted by |T|p_p. We will need to compare this inner product, 
with its metric g, to one using the Euclidean metric, which we denote by S, where S^'^ is the 
Kronecker S. Thus, for the Euclidean case, which depends on the chart (il, ^), we have 

i,iG{l...d}'= ie{l...d}>' 

with |T|5^2. being the corresponding norm. In the sequel, we will need the following lemma com- 
paring lengths of tensors in these two inner products. 

Lemma 2.1. Let x & U = ^(it). If Ax{g) and Xx{g) are the largest and smallest eigenvalues of 
g^'^{x), then 

Xa,{gf{T, T)5,. < (T, T)g,p < A,(5)'(T, T)^,,. (2.3) 
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Proof. The matrix g^'^ is positive definite, and so each eigenvalue Xm is positive. If necessary, we 
repeat an eigenvalue. Let be the corresponding eigenvector. The set of eigenvectors is assumed 
to be orthonormal. If it^ is the i*'* entry in the column u^, the spectral theorem implies that 
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m£{l...d} m£{l...d} 

Using this in connection with the expression for (T, S)g^p yields 

/ V 



(T,T),,p 



^ ^ Xmi ■ ■ ■ A 

) 

/ 



rufe 



(mi...mfc) 
G{1...4'= 



\e{i...rf}'= 



^ Aa;(5() n^^tt:^^ ■ ■ ■ n^^ti:^^Tj^^...^i^Tj^^...jj, 

(«i---«fe) (ji-jk) {mi...mk) 
Ve{l...d}''' G{l...rf}'= GiL-.d}* / 

A.(5)' E (Tn,....j' = A.(5)'^(T,T)5,., 

G{l...d}'= 



which establishes the lower bound in (j2.3p . The upper bound follows in the same way. 



□ 



Geodesies A curve 7(t), defined by (x-'(t)) in local coordinates, has its tangent vector given by 
7 = X^jgji d} ^^j- '^^^ length of 7(t) is defined by 



If p,q are points in M, then the distance between p and q, dist(p, g), is defined to be the infimum 
of the length of all piecewise differentiable curves joining p and q [3 §7.2]. 

Extremals of the arc length functional are called geodesies. If we use the arc length s as the 
parameter (i.e., t — )• s), then the Euler-Lagrange equations for x^{s) are the following second order, 
nonlinear differential equations: 



d'^x'' -sr^ pfc dx"^ dx^ 



0. 



(2.4) 



The are called Christoffel symbols; they are given by 



2 E 

me{l...d} 



km 



dgjm , dgim dgi 



+ 



dx^ dxi dx^ 



(2.5) 
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A geodesic solving (j2.4p is specified by giving an initial point p G M, whose coordinates we may 
take to be (x^(0),... ,x^{0)) = 0, together with a tangent vector tp having components ^(0). 
A Riemannian manifold is said to be complete if the geodesies are defined for all values of the 
parameter s. All compact Riemannian manifolds are complete. M'^ is complete. 

Exponential map We define the exponential map Expp : TpM — > M by letting Expp(O) = p and 
Expp(stp) = 7p(s), where 7p(s) is the unique geodesic that passes through p for s = and has a 
tangent vector 7p(0) = tp of length 1; i.e., (tp,tp)gp. = 1. By the Hopf-Rinow Theorem [9l §7.2] 
M is complete if and only if Expp is defined on all of TpM. Other important consequences of this 
theorem are that M is a complete metric space under the distance dist(p, g), that any two points 
p,q €M. are connected via a geodesic of minimum length dist(p, q), and that the exponential map 
is defined on the whole tangent bundle TM; as such, it is a smooth function of both arguments 
la §9.3]. 

Although geodesies having different initial, non-parallel unit tangent vectors tp = 7p(0) may even- 
tually intersect, there will always be a neighborhood Up of p where they do not. In ilp, the initial 
direction tp, together with the arc length s, uniquely specify a point q via q = 7p(s), and the expo- 
nential map Expp is a diffeomorphism between the corresponding neighborhoods of in TpM and p 
in M. In particular, there will be a largest ball i?(0, rp) € TpM about the origin in TpM such that 
Expp : B{0,rp) h{p,rp) C M is injective and thus a diffeomorphism; rp is called injectivity radius 
for p. By choosing cartesian coordinates on i?(0,rp), with origin 0, and using the exponential map, 
we can parametrize M in a neighborhood of p via q = Expp(x), x € TpM. Finally, any chart with 
11 = expp(i?(0, r), r < rp and 4> = Exp"^ is called a normal chart of radius r and corresponding 
cartesian coordinates are called a normal coordinates. 

The injectivity radius of M is im '■= infpeAf ^p- If < rM < oo, the manifold is said to have positive 
injectivity radius. All compact Riemannian manifolds without boundary are both complete and 
have positive injectivity radius. On the other hand, many noncompact Riemannian manifolds, 
including the simple case of M"^, also have positive injectivity radii. In fact, the injectivity radius 
of M'^ is oo. It is worth pointing out that manifolds having rM > are always complete [13^ §1.1] 

We make special note of the fact that, for a compact Riemannian manifold, the family of exponential 
maps are uniformly isomorphic; i.e., there are constants < Fi < r2 < cxd so that for every pQ £Wl 
and every x,y £ B{0, r) 

Ti\x -y\< dist(ExppQ(x),Expp,,(?/)) < T2\x - y\. (2.6) 
We will have more to say about this later in Remark 12. 4[ 

2.2 Covariant derivative bounds 

We now turn to the concept of higher order covariant differentiation on manifolds. Such derivatives 
are intrinsic; they are used to define Sobolev spaces on measurable subsets of M. After relating 
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covariant derivatives and ordinary partial derivatives in equations (j2.7p and (|2.9p , we obtain bounds 
that will be used to prove Lemma 13.21 which is a key result concerning local norm equivalence 
between Sobolev spaces on the manifold and Sobolev spaces on normal charts. 

The covariant derivative, or connection, V associated with (M,(7) is defined as follows. Let v = 
^"j^i v^^j and w = X^j=i w^ej be vector fields. The covariant derivative of w in the direction of 
the vector field v is given in terms of the Christoffel symbols, 

Vvw:= E E [(.^^+ r^-W)]e„ 
fce{i...d}iG{i...d} 'je{i...d} 

and is itself a vector field. If we use a covariant vector field, w = Ylke{i d) ^fc^'^i then 

Vvw= E E [K^- E r>V)]e^ 

feg{l...d} je{l...d} rG{l...d} 

This is a directional derivative. The appropriate "gradient" is a rank 2 tensor, 

(fc,j)e{i...d}2 \ se{i...d} I 

Covariant derivatives may be defined for any type of tensor field - covariant, contravariant, or 
mixed. For instance, if T is the order k (covariant) tensor defined by 

T= ^ Tje*! «) • • • ® 

i&{\...dY 

then the covariant derivative of T is 

(dT- ^ 
r=l se{l...d} 

which is an order k + \ covariant tensor. Of course, if T = / is a scalar valued function, then 
V/ = Eje{i...d} (oi'' ^o'^e precisely, V/ = Ejg{i...d} e^) which is just the usual 

expression for the gradient. The "Hessian" is V^/; it has the form 

^ ^ dx^dxo ^ ^'-^ dx^ 

(i,j)G{i...d}2 V fce{i...4 

The components of the A;*^ covariant derivative of / have the form 

fc-i 

(vV(x))j = {d^!{x)\ + E E (2-7) 

m=l Jg{l...d}™ 
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where 

and where the coefficients x i— )• A-l{x) depend on the Christoffel symbols and their derivatives to 
order k — 1, and, hence, are smooth in x. This can also be written in standard multi-index notation. 
Let ai, a2, ■ ■ ■ , ctd be the number of repetitions of 1, 2, . . . , d in j, and let a = (ai, . . . , a^). Then, 

^^"^^^^ ■= dix^r^ • • • dix^^ ^ ° =■ ° ^''^ \a\ = Y:»k = m. (2.8) 



Unlike ordinary partial derivatives {d^f{x))i or even the Christoffel symbols, the components 
(V^/(x))j (and, of course, those of VT) transform tensorially under a change of coordinates, 
so that V^/ is independent of the coordinate system. This was a discovery of Levi-Civita [9, §2.3]; 
it is why the covariant derivative V is an important improvement on the partial derivative d. 

Even though the partials are not invariant, we can express them in terms of covariant derivative 
components, provided we use local coordinates. Because of the relationship between the highest 
order derivative in (j2.7p and the corresponding covariant derivative's component, we can solve for 
the highest derivative in terms of covariant derivatives: 

{d'fix))i = {V'f{x)),+ J2 E Bi{x)iV^f{x)),„...,^. (2.9) 

m=l 



The expressions in ()2.7p and ()2.9p are needed to obtain bounds on (V'^/, f)g^p. These bounds 

will prove useful in estimating certain invariant Sobolev norms using local coordinates. They are 
given below. 

Proposition 2.2. Let (il, </>) he a chart, with U = 4>{ii). For every x = (f){p) G U, there exist 
multivariate polynomials cu{g,k) > 1 and Cu{g,k) > 1 in gij{x), g^^ {x) and their derivatives up to 
order k such that 

k k 

IVVIL < Cu{g, k) \d'^f\l. and \d'f\l, < cu{g, k) ^ | V-/||,, (2.10) 

m=l m=l 

where \ ■ \s^x is the norm associated with the inner product in 12. In addition, we have that 

k k k 

cu{g, k) J2 Id^'flh < E l^'/I^.P ^ ^^(5, k) J2 I^'^/IL^ (2-11) 

m=l m=l m=l 

where cu{g,k) = {Y^'^^^mm{l, Xx{g))-"^cu{g,m)))~^ and Cu{g,k) = Ylm=i^x{g)"^Cu{g,rn)); 
Xx{g) and Ax{g) are the smallest and largest eigenvalues of g^^(x). 
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Proof. From (j2.7p and Schwarz's inequality, we see that 

|(V'=/),p < A? (\{d'f f + J2 \d"'f\h] 

\ m=l / 

Where ^? := 1 + E™=i E^e Summing over the i indices yields 

k 

IVVIL < Cuig, k) ^) := E A- (2.12) 

Examining the coefficients in (j2.7p shows that Cu{g,k) is a polynomial in the Christoffel symbols 
and their derivatives. But these are linear combinations of derivatives and multiples of the metric 
components, so Cu{g, k) is a polynomial in g^^ , gij and their x derivatives. Repeating the procedure 
with (j2.9p gives the result for \d^f\'j^- The inequalities in (|2.14p follow from (|2.13p and Lemma [2.11 

□ 



Let q G M. Suppose < r < rM and sufficiently small that, in a fixed chart (11, (/>), il contains 
a neighborhood of q and all of the associated geodesic balls b(g, r). As we mentioned earlier, the 
exponential map is a smooth on the tangent bundle. Consequently, if we let x,y = 4>{q), € U , then 
p{x,y) = Exp^-i(j^)(x) € M is smooth in both x and y. If q is allowed to vary, then in a normal 
chart (b(g, r), Exp~^) the entries of the metric tensor gi j, its maximum and minimum eigenvalues, 
its determinant, and the Christoffel symbols vary smoothly in both y = <j){q) as well as x. They 
are thus uniformly continuous on the closure of Ugb(g,r), and are bounded there independently of 
q and x. Since we can cover the compact manifold M with a finite number of such neighborhoods, 
the boundedness of the quantities associated with the metric tensor hold uniformly for any normal 
coordinate system of radius r. 

The manifold M doesn't have to be compact for the metric and related quantities, when expressed 
in normal coordinates on b(g, r), to be bounded uniformly in q by constants depending only on r. 
This also holds for noncompact C°° Riemannian manifolds having hounded geometry^ which means 
that they have positive injectivity radius and that for all fc > 0, the curvature tensor R satisfies 
suPpgM \'^^R\g,p < Cfc) where Ck depends only on /c [U [12l [131 EE]- Compact manifolds fall into 
this class, as does the space W^. The implication that is important here is the following result. 

Corollary 2.3. Let M 6e a compact manifold, or, more generally, have hounded geometry. In 
addition, suppose < r < rM and H = b(g,r) and (f> = Exp~^. Then there are positive constants 
c{i,k), C{i,k), c{i,k), and C{i,k) such the following hounds hold on Exp^^(b(q, r)) independently 
of q: 

k k 

\V'f\l, < (7(r, k) E \d"^f\l, and < c(r, k) \V"^f\l,, (2.13) 

m=l m=l 

where \ ■ \s^x is the norm associated with the inner product in \2.S^) . and 

k k k 

c(r, k) \d'^f\h < E l^'/l'-P ^ ^(r' ^) E I^'/IL- (2.14) 

m=l m=l m=l 



13 



Proof. All of the constants in cu{g, ni)) and Cjj{g, k) in (j2.10p are polynomials in various quantities 
associated with gij. Since these quantities are all uniformly bounded by functions of r, we see that 
there are new constants c(r, h) and C(r, k) such that cu{g, k) < c(r, k) and Cu{g, k) < C(r, k). The 
bounds in (]2.13p follow similarly. □ 

Remark 2.4. From what we said earlier, we note that in a normal chart of radius r < rM, the 
measure on M has the form = ^ det {gij)dx^ ■ ■ ■ dx'^. If M is compact or has bounded geometry, 
then there are positive constants ci(r),C2(r) such that ci(r) < det (gij) < C2(r) on i3(0,r). We 
also point out that, for similar reasons, Fi and r2 in ()2.6p depend only on r, and not the center po. 



3 Sobolev spaces on subsets of M 



Sobolev spaces on subsets of a Riemannian manifold can be defined in an invariant way, using 
covariant derivatives [2] . In defining them, we will need to make use of the spaces Lp, Lq. To avoid 
problems with notation, we will use the sans-serif letters p, q, rather than p, q, as subscripts. Here 
is the definition. 

Definition 3.1. Let Q C M. be a measurable subset. For all 1 < p < oo, we define the Sobolev 
space Wp''{il.) to be a// / : M — )• M such that |V/|g^p in Lp. The associated norms are as follows: 

, (E^=o^lvVl.%dMp))'/^ i<p<oo; 



maxo<fe<m|| I V^/lg^pll^^j.^^, p 



oo. 



When p = 2, the norm comes from the Sobolev inner product 

m „ 

{f.g)m,n := {f,g)w,-(n) E / (VV, V^ff) d/.(p). (3.2) 

We also write the p = 2 Sobolev norm as WfW^Q ■= {f,f)m,n- When Q = M, we may suppress the 
domain: (/, g)m = (/, g)m,M and 



When is a region in Euclidean space M"^, the definitions above are equivalent to the standard 
ones (cf. [H Chapter III]). In particular, when p = 2, the definitions coincide. Furthermore, if the 
p = 2, and the M.'^ norms are expressed in the multi-index notation in (12. Sp . then they have the 
form 




3.1 Metric equivalences 



There is a metric equivalence between the M'^-Sobolev norm on C -6(0, r), where r < rM and the 
M-Sobolev norm on the image Exp (il) C h{pQ,r), pQ G M. In particular, the following lemma 
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shows that the exponential maps induce operators / i— )- / o Exp^^ that are boundedly invertible 
from all spaces l^™(Expp(J7)) to W^{Q). This is especially important for us because the constants 
involved are dependent only on r, and not on po; that is, the equivalence between these spaces is 
independent of pq. 

Lemma 3.2. For m E N and < r < i'm, there are constants < ci < C2 so that for any 
measurable Q C Bj-, for all j € N, j < m, and for any po £ M, the equivalence 

cilluoExpp^ llty.(n) < \Mw^(E.p^^{n)) ^ C2||noExpp^, il^,^^^ 

holds for all u : Expp^{il,) M. The constants ci and C2 depend on r, m and p, but they are 
independent of and pQ. 

Proof. For Q C B{0, r) and 1 < p < oo we have 

II^III Exp,^,(n) = E / |VVlg%V^et(5,,(x))dx. 



By Remark 12.41 yjdet[gij{x)) is bounded above and below by constants depending only on r. 

Similarly, by Corollary [231 the pointwise metric norms of the covariant derivatives \^^f\g,p are 
bounded uniformly by the corresponding Euclidean quantities. Combining these facts yields the 
result for p < cx). The p = oo follows similarly. □ 

3.2 Sobolev bounds for functions with scattered zeros in M 

In this section, we show that Sobolev space functions having many zeros are uniformly small, pro- 
vided the Sobolev space has high enough order. An immediate consequence of this gives pointwise 
error estimates for many kinds of interpolation processes. We will discuss this later, in Corollarv l3.91 

The Euclidean case We will need to discuss bounds on certain Sobolev norms for a special 
class of domains in M.'^. Following Brenner and Scott [U Chapter 4], we will say that a domain 
D is star shaped with respect to a ball B{xc,r) C D if for every x G T>, the closed convex hull of 
{x} U B{xc, r) is contained in 2). If T) is bounded, then there will be a ball B{xc, R) that contains 
D. Of course, the diameter dx) of T> satisfies 2r < d-D < 2R. An important, useful geometric 
quantity associated with D is the chunkiness parameter 7, which Brenner and Scott [3^, Definition 
4.2.16] define to be the ratio of d'Q to the radius of the largest ball -Bmax relative to which T) is star 
shaped; i.e., 7 = du/rmax, where rmax is the radius of -Bmax- It is easy to see that 7 < In the 
case T) = B(xc,r), which is star shaped with respect to itself, the chunkiness parameter is 7 = 2. 
Another geometric property associated with D is a cone condition. Every x in 2? is the vertex of 
a cone Cx € D whose axis is along x^. — x, radius is r, and aperture (maximum angle across the 
cone) is 9 = 2 arcsin (275) [25, Proposition 2.1]. 
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Lemma 3.3. Let X = {xi, . . . ,xiy} C D, lefN 3 k > 0, and let h > 0. Suppose that h = hx, the 
mesh norm of X in D, or, less stringently, that every hall B{x,h) C D contains at least one point 
in X. If h satisfies 

rsm{9) 

and if p G 7rfc(R'^), then for every multi-index a, \a\ < k, then 

2k^ 



\D ^'llioo(D) < 2 ( — ^ ) \\p\\i^(x)- 



Proof. Apply [25", Proposition 2.3], as modified by [251 Remark 2.4]. □ 

Remark 3.4. Two things. First, if /i = hx, tfien every ball B{x,hx) C D contains a point in X. 
Thus the "or" is not necessary. Second, if we take r = rmax and R = dri, then ^ = = 27 and 

' 'max 

sin(^/2) = It is easy to show that 

rfD rsin(g) 2fc^ 47^fc2 

- 4(1 + sin(e))fc2 rsin(e) " dj) ' ^ ' 

Thus the restriction (j3.3p on /i and the bound on |1^°p||l^(d) now become 



Proposition 3.5. Let D dW^ he a hounded, star-shaped domain, m € N and p G M, 1 < p < 00. 

Assume m > d/p w/ien p > 1, and m > d, for p = 1. If u ^ W^{T)) satisfies u\x = 0, where 



X = {xi, . . . , xat} C 2) and if h = hx < .jjS-i , then 



W\w^{D) ^ Cm,d,pl dl^ \u\w^{v)- (3-6) 
II II ^ /-< djni~ii/P| I /o ^\ 



Proof. Following the notation in Brenner-Scott [3], we take Q"^u € TTm~i{^'^) to be the Taylor 
polynomial for u averaged over i?max) the largest ball relative to which 2) is star shaped. We begin 
by estimating \Q"^u\yyk^^-j))- Using u\x = 0, we have 

By [31 Proposition 4.3.2], with the 7-dependence involved there explicitly included, the right-hand 
side above has the bound 

WQ'^u - u\\l^(^x>) < Cm,d,pl'^d'^~'^^^\u\w;^('r,y (3.8) 

It follows that 
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The next step requires applying Lemma [33l with r = rmax and R = cId, and then using Remark [3.4[ 
Doing so results in 

iin"r)™o,ii <r r<' ^,'^+2|a|^'^-l"l-'^/Pu,i fio^ 
11-^ Q ^IIloo(D) < '-'m,d,p7 '"b \u\w^iv)- W-yj 

Sum over \a\ = k < m and use ||L'"Q™n||^^fe(j,) < ^^/''[[^"^'"'"[[^^^(d) to get 

|/^m I ^ /-(// d+2k jm—k\ i /o i n^ 

We need to estimate \u — Q"^u\y^k(j)y To do this, we will use the Bramble-Hilbert Lemma [3], 
which holds for 1 < p < oo. Take r = rmax in the lemma. We have, after carefully tracking the 
7-dependence of the constant there, 

\u-Q"'u\wk^^)<Cm,d-f'^d:^-''\u\w^m^j:)), k = 0,1, . . . ,m. (3.11) 

Using the triangle inequality in conjunction with (|3.1U|) and (j3.1ip results in (|3.6|) . In addition, 
doing the same with the bounds in (j3.9p . for |a|, and ()3.8p yields ()3.7p . □ 



The inequality (|3.7p is a special case of the one established next. 
Corollary 3.6. With the notation of Proposition [375[ we have 

Proof. From our estimate on |(5'"'u|vi/4(D)) it is easy to show that 

Putting the two estimates together via the triangle inequality yields p.6p . Now suppose that 
u G W^(a?). D^Q^n = so 

||Z)"(Q™n-^.)|| = IIQ-^I-Ip^^x) - (p'-u)\\ < C„_|,|,,,p(l + 7)'dr'"'"'^'KI^;n-|.|(25)- 
Hence, we have 

l/nm I ^ r< djm—k—d/pi i 

^ ~ "^^114^4(1)) — '-^m-fc,d,p7 "d l^lH/p'"^''(D) ' 

By (|3.6p . k replaced by m — A;, we finally arrive at 

Combining this with our bound on \Q^u\y^k and again employing the triangle inequality, we 
obtain the desired inequality. □ 

We will next apply the result above in the special case of a ball, where the bounds simplify con- 
siderably. Specifically, when u G W™(i?(a;,r)) and u\x = 0, it allows us to control certain sums of 
lower order Sobolev norms by |w|w''"(B(a;,r)) • Doing this yields the following: 
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Lemma 3.7. Let X = {xi, . . . ,X]\i} C B(x,r) have its mesh norm h = h{X,B{x,r)) satisfy 
h < hor, where ho := In addition, suppose that m > d/p, if p > 1, and that m > d, if p = 1. 

Then there is a constant Cm.d.p > such that the estimate 

/ ^ rP(^"'"Vl^fc(B(x,r)) ) <Cm,d,pHw^{B{x,r)) (3-12) 

holds for all u € W^{B{x^r)) vanishing on X (i.e., = Oj. In addition, we have that 

\W\\Lac{B{x,r)) < CTn,d,pr'^~'''^'^\u\w^{B{x,r)) (3.13) 



Proof. For a ball D = B(x,r), which is of course star shaped, the chunkiness parameter is 7 = 2, 
the diameter dB(x,r) = 2r. Since h < ^Q^^^i = hor, Proposition 13.51 applies. Thus, for A; = 0, . . . , m, 
the bounds in ()3.6p become 

Standard algebraic manipulations of the expression above then yield p.l2p . The last inequality 
(j3.13p is a direct consequence of (j3.7p . □ 



The manifold case The case of Sobolev bounds on u when the underlying set is a geodesic ball 
in M can be treated using a combination of the results involving metric equivalence. Lemma 13.2^ 
and the corresponding Sobolev bounds in Lemma [3. 71 for Euclidean balls. We can treat much more 
general situations than the one described below, but for now it is precisely what we need for the 
sequel. 

Lemma 3.8 (Zeros Lemma). Let m he a positive integer, greater than d/2, and let r he a positive 
real numher less than tm, ihe injectivity radius o/M. Suppose that H C h{p,r) C M. is a discrete 
set with mesh norm h < Lir/iQ. If u G W|"(b(p, r)) satisfies u^^ = 0, then for every q G h{p,r), 

\u{q)\ < Cm,Mr'""''/^||u||vF2'"(b(p,r)): 

where Cm,M is a constant independent of u,p,q,h and r, and where Fi is as in ^2. 



Proof. Using the diffeomorphism Expp, we set X = Exp~ (H) and note that this set has mesh 
norm h{X,B{{),r)) < rho by (12. 6[) . Defining S as u o Expp, and setting q = Expp(2:), we see that 
(|3.13p applies to u, giving 

u{q) = u{z) < Cm,dr^"'''''^'^\\u\\w^(^B{0,r)) < Cm,M''^'^~'^^'^\\''^\\w^(b{p,r))^ 

by Lemma |3.2[ □ 
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3.3 The family of kernels Km,M on M 

We now are prepared to identify the family of kernels associated with bounded Lebesgue constants. 
A well known fact, which also happens to be a simple consequence of combining Lemma 13.21 and 
the Sobolev embedding theorem on domains in M"^ via, is that W|"(M) is embedded in the space 
of continuous functions on M, for m > d/2 [21 §2.7]. Consequently, point evaluation is a bounded 
linear functional, and has a unique reproducing kernel, which we define by Km,M : M x M — >■ M, 
although we often suppress the domain, writing Km = ^m.M- Being a reproducing kernel means 
that 

for all / € W^. The reproducing kernel is necessarily strictly positive definite: the formula 
Cgh ^?^C^™(^' ^) = II X^^gs C)llm = implies that there exist coefficients {a()(£'E so 

that, for all / G W^, /(^) = l^fgH\{^} '-'^c/(0- Using a bump function centered at ^ for / easily 
provides a counterexample. 

As an aside, we note that we can modify the Sobolev norms (and, hence, the reproducing kernels) 
in the following benign way: 

m „ 

A:=0 

where Cm > 0, Co > 0, and Ck > 0, k = 1, . . . ,m — 1. For such a modified norm. Lemma [312] holds 
in precisely the same way, except with different constants ci, C2. The benefit, when M = M"^ is that 
for a particular choice of constants Ck we have the inner products corresponding to the Sobolev 
(or Matern) splines [22j . 

K{x,a) = Cm,d\xr''^^^Kd/2-m{\x - a\) 

where is a modified Bessel function. This is achieved when ||/||^ ^^d^ = f^d fix){l — ^)"^f{x)dx, 
which is easily accomplished because Jj^d /(x)A™/(x)dx = (— 1)™ /jgd(V™/(x), V'"/(x))dx holds 
by integration by parts. 

Positive definite kernels and their native spaces The situation above can be turned around, 
in the sense that we are able to start with a symmetric, strictly positive definite kernel k, and 
then construct a corresponding Hilbert space that is the RKHS for it. The term strictly positive 
definite means that for any finite set H C M the interpolation matrix 6= = (^(Cj 0) ^)g=2 is 
positive definite. The construction of the corresponding RKHS for such a kernel is described in 
detail in [29, §10.2]. The space itself is known as the native space 3sf(K), and its inner product is 
denoted by Because C= = (k(C, 0) ^)g=2 is positive definite, the interpolation problem 

S|„ = /|_, for s G span^g2^(''C) always possesses a unique solution, denoted by /=/. Equivalently, 
/=/ can be determined by finding the solution to the variational problem argmin{|[s||>j-(K) : s G 
?vf(K), S[„ = The relationship between positive definite kernels and their native space inner 
products translates into a duality between kernel interpolation and the variational problems with 
interpolatory constraints. This is discussed in detail in the celebrated Golomb- Weinberger paper 

m- 
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For a kernel associated with a radial basis function, the native space has its origin in the work 
of Madych and Nelson (see [2T] for an example), where the kernel interpolation problem is recast 
as a variational problem, one where the interpolant is the minimizer of a Hilbert space norm (or 
seminorm) over all possible interpolants of data. The native space appellation itself is due to 
Schaback [27] . who extended this idea to treat kernels on more general domains. 

We now turn to native space error estimates for interpolation by positive definite kernels associated 
with Sobolev native spaces. The following corollary to Lemma 13.81 shows that, regardless of the 
kernel giving rise to a Sobolev native space, it is always possible to measure the approximation 
order for interpolation from the native space. This generalizes the previous native space result for 
R'^ and S<^. 

Corollary 3.9. Let k : M x M ^ M : (??, C) ^ ^^(^1 C) ^6 a positive definite kernel with native 
space ?sf(K;) = VF™(M), m > d/2. Assume that H C M has mesh norm h < iMhoTi. Then for 
f G WY^(M.), the error incurred by interpolating with k at the nodes H C M is 

\f{x)-Hf{x)\ < CMh'^-''/^\\f\\w.rm- 

Proof. By picking r = h/{hoTi), we have r < rM and Lemma 13.81 applies to the ball b(p, r), giving 

Since, by assumption || • |Ivi/2™(m) ^ II ' lU) we have the chain of inequalities ||/ — Isf\\w^(M) ^ 
C||/ — /s/IU ^ C'II/IIj^ < C'||/||iy2™{M); where the middle inequality follows by the Pythagorean 
theorem, since is an orthogonal projector on the native space. □ 



4 The Lagrange Function 

In this section, the Lagrange function centered at an arbitrary point ^ G S (usually suppressing the 
subscript X = Xs.) is investigated. We begin by showing that for a fairly general class of kernels, 
X is bounded for quasiuniform centers. In Section [4.21 it is shown shown that, for a specific class of 
kernels, x is actually controlled by a rapidly decaying function of dist(x, Specifically, in the 

ball b(^,rM) about ^, \x{x)\ is controlled by exp[— i/dist(x, (this is Proposition 14. 5p . 

4.1 Uniformly Bounded Lagrange Functions 

Our first goal is to obtain bounds on the decay of the Lagrange (or fundamental) function for kernel 
based interpolation. As before, we denote the native space for a positive kernel by [N"(k). 

Definition 4.1. Given a positive definite kernel k : — )• M, and a finite set E C M. we denote the 
Lagrange function centered at^ €E by xs,- Le., X(,{C) = ^i^-, C) fori, C ^ ^ o-^idxE^ ^ span^g^ ^('j 0- 
By the discussion in Section 3, we see that Xi = argmin{||s||jj-(K) : s{C,) = S{i,C,),C, G H}. 
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We first observe that, when the centers are quasiuniform, the function x is bounded. 

Lemma 4.2. Suppose that k is a positive definite kernel on M with native space = W^(M), 
m > d/2. If the centers H are quasiuniform - namely, there exists a constant p such that h/q < p, 
- then the Lagrange function x is bounded, with a constant depending on k, m and p only. 

Proof. Let ijj he, a, 'bump' function with support inside the ball b(^, g) = {a € M : a) < q} 
obtained by dilating a univariate function: 

where a € C°°([0,oo)) is decreasing, has support in [0,1), and satisfies (t(0) = 1. The Lagrange 
function interpolates ip, and Corollarv 13.91 provides the estimate: 

||V'-xlloo<c/i™-'/'||^IU. 

From its definition, it is clear that ip is bounded: \\ip\\oo is 1. Consequently, 

||xl|oo<l + C/l"-'^/2||^|U. 

Thus, we need only to estimate this last quantity. Since ^ is a function of the distance from 
^ only, we can use Lemma 13.21 to estimate it. Note that, if Expg(x) = a then ip o Exp^(x) = 
^ j- dist(Exp,(.).Exp,(o)) N^ = (t)- So Mm is Controlled by \\a (^) ||^ < which means 

that 

llxlloo < l + C/>"^-^/2. 

Here C is the constant determined by the interpolation error, the constant from Lemma 13.21 and 
the m^^ Sobolev norm of the univariate function a. □ 



4.2 Rapid Decay of Lagrange Functions 

For the kernels K,m,M, we can improve significantly over Lemma 14.21 by employing an argument of 
Matveev [23], to obtain very rapid decay of the Lagrange function. The strong metric isomorphisms 
provided by (j2.6p and Lemma [3^2] permit us, roughly, to ignore the manifold and carry out most of 
our analysis on the tangent plane. 

This has the additional benefit that, although it is not compact, one may take M = M*^. In this case, 
the Riemannian distance is simply the Euclidean distance dist{x,y) = \x — y\ and the exponential 
map at is nothing more than translation by ^: Fjxp^{x) = x + ^. Furthermore, the injectivity 
radius is rM = oo, the constants from (12. 6p are Fi = r2 = 1, and, likewise, the constants from 
Lemma 13.21 are ci = C2 = 1. The underlying kernel is the one associated with the Sobolev norm on 
M!^ and is very similar to the so-called Sobolev (or Matern) kernels. In fact, a benign modification 
to these arguments, namely, by replacing the Sobolev norm with either the (previously defined) 
Sobolev seminorm (in fact, this is the variational problem originally considered by Matveev [23]) 

\f\l,n= [ (V-/(x),V"/(x))dx 
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or with a slightly reweighted Sobolev norm 



2 

k=0 



m „ 
7 — n "'^ 



as discussed in Section 13.31 gives totally equivalent results for the surface splines and the Sobolev 
splines, respectively. 

The pointwise Lagrange function estimate is essentially a "bulk chasing" argument, adapted from 
|23t Lemma 5], where the central observation is that the "bulk" of the tail of the Lagrange function's 
Sobolev norm is always contained in a narrow annulus. This is a property shared by exponentially 
decaying functions: the integral of the tail of a nonnegative function of this type can be controlled 
by the integral on a sufficiently wide annulus. The next lemma makes this precise: 

Lemma 4.3. Assume that S C M has mesh norm h < Fi/iq min(l, rM/3). There exists a constant 
€ E (0, 1) (depending on d and m) such that for 1 <t < ''"3^^^ 

ll^li„..bC(,,3<^)) ^ ^II^IL,b^(,,3(*-l)(^))' (4-1) 

where Ti is the constant from the left hand side of \2. 0\) . xm is the radius of injectivity and ho is 
the constant from Lemma \^.7\ 



Proof. We start by constructing a C°° cutoff, (j), that vanishes outside b(^, (3t — and which 

equals 1 on b(^, {3t - 2)p^). Let 

(j)(a) = a ( ^4^dist(a,0 - 3(t - 1) 



where (t(T) equals 1 for T < 1 and vanishes for T > 2. Since (p is a function of the distance from ^ 
only, we can rewrite it as a composition of a univariate function and Exp^. Let Expg(a;) = a, then 

(^oExp5(x) = C7 (^^^dist(Expg(x),Exp^(0)) -3(t- 1)^ = a (^^\x\ - 3{t - 1)^ . (4.2) 



The region where (f) is nonconstant - a(^,3t — l,p^) - is the middle third of the annulus 
a(^,t,3p^). In order to simplify notation, we'll denote the annulus where (p transitions from 
1 to by the gothic letter a. That is. 

Likewise, we'll denote the region where (p equals 1 by 
and by b'' its complement M\b. 
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By the variational property of X) the Sobolev norm of x is less than that of </>x, since the two 
functions coincide on H. Moreover, because of the compact support of 0, the Sobolev norm of (px 
can be expressed in terms of its behavior on b U a, noting that on b it is identical to x- Thus, 
we have that Hxllm — llxllmfa + ll<^xllmo- -^y subtracting and applying Lemma [321 we obtain the 
estimate 

where we identify corresponding annulus of interest in T^M by 

>l = ^(0,3t-l,-^). 

i itlQ 

The annulus A is the preimage, via the exponential map, oi a, x'-= X° Exp^, and cp = (p o Exp^. 




Figure 1: Important annuli on the tangent space T^M. 



Our goal for the rest of the proof is to estimate the quantity ||0xllmyi- This involves two key 
estimates. The first estimate is concerned with the effect of multiplying by the cutoff (which we 
can tackle by using Leibniz' rule). The effect is independent of the radius t, and allows us to 
control this norm in terms of a combination of seminorms of x- The second estimate will control 
this combination of seminorms via partitioning the annulus into small caps where Lemma 13.71 can 
be applied. 

Estimate 1 The desired result is to control 110 xWmA ^ biased Sobolev norm of x on yi - this 
is g31). 
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We observe, first of all, that we can simplify matters by applying ()4.2p and the product rule: 

2 



'Xllm,yl 



E 

7I 



1 J Ja 



/3<7 



dx 



|7|<m /3<7 



^,IjA)|^l_3(t-r 

h 



dx 



A direct application of the inequality, |L''''(T(r|x| — T)| < Crl'^'l , which holds for T > 0, 
shows that is bounded by Cm,d^\^\<m IaT^^k-j (J^) ^ \D^xix)f dx. Because 

Y~f^ < 1, and f " j < f \° j , we are left, after rearranging terms, with the estimate 



< 



Cm,d ^ ^ 



|/3|<m 



r,/,o\2(— 1/31) 



A 



dx 



(4.4) 



Estimate 2 We are now in a position to apply Lemma 13.71 We cover A with a sequence of balls 
[Bj)j^^ such that 

• each ball is of radius 



each ball is contained in the large annulus ^(0, f,3p^) (displayed in Figure 1). 

every x E j4(0,t, Sjr^) is in at most N{d) balls Bj, with N{d) depending only on the spatial 
dimension and not on t, h, hQ,d, etc. 



Observe that the combination of seminorms from (j4.4p can be bounded by corresponding norms 
carried by the balls Bj. The mesh norm in each Bj is h/Ti < ho (at most), so Lemma 13.71 applies. 



with r 



E 

|/3|<m 



^ \2{|/3|-m) 



Li/io 



Li/i, 



2{|/3|-m) 



< Cm.M E I ^ I H/j'" {-Bj ) 
■|2 



< iV(o!)C„,M|xl 



iyr(A{0A3p^)) 



The second inequality is Lemma 13.71 while the third inequality is a consequence of the finite 
intersection property of the balls Bj . By a direct application of Lemma 13.21 we obtain 



|/3|<m 



h 



2(l/3hm) 



\D 



A 
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Putting this together with (j4.4p and (j4.3p . we see that the Sobolev norm of the Lagrange function 
taken over the complement of a ball can be controlled by the Sobolev norm on a thin annulus near 
its boundary: 

< Cm,M\\xf \=-K\\xf ft V (4.5) 




Figure 2: Important sets on the manifold. 



What remains, is to reinterpret this as an inequality involving only complements of concentric balls. 
It follows from their definitions that b = b (^^, (3t - 2)p^^ C b (?,3tp^^ =: b+. This implies 

that (b^)^ C and, setting b^ := b 3{t — l)p^^ , we make the simple but useful observation 
that a(^,t,3p^) = {b~)^ \ (b+)^. Applying this to ([45]) we get 



and, consequently: 



I ||2 < K I W 11^ _ II ||2 

llxlL,bC(^,3ij^) ^ ^II^IL,bC(c,3{t-l)i^) 



with e := a/— ^— < 1. □ 
l + K 



For a sufficiently large value of t less than a fixed multiple of the radius of injectivity, the previ- 
ous lemma could be repeated several times, to obtain better estimates, resulting in the following 
Corollary. If we have M = R"^, the radius of injectivity is rM = oo, and t may become arbitrarily 
large. Furthermore, in the Euclidean case one may choose the surface splines as the kernel (instead 
of KM,m)) which results in [23^ Corollary p. 130]. Alternatively, one may choose the Sobolev splines 
(see Section 13. 3p and achieve precisely the following: 
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Corollary 4.4. There is a constant > such that, if 'E has mesh norm h < Fi/iq min(l, im/S) 
then for < T < im we have the estimate 



\\X\\ni,hH^,T) ^ '-m,Me ' h ' \\X 

where q is the minimal separation distance between points of 



d/2-m 



Proof. The corollary is a consequence of iterating Lemma [4. 3 1 as often as possible. Let T = Stp^]^ = 
3np^ + r with < r < Sj^- The inte ger n represents the number of times the inequality ()4.ip 



can be iterated. It follows that 



\x\L,hH^,T) ^ e||xlL,bC(5,3(t-i)T^) 



< e"||xl|m < e ^e* ||x||m 
_i / Eihsi\ii II II 

The first inequality in the statement of the corollary follows with = |^-^ loge| and is absorbed 
into the constant. 

The norm of x can be estimated by using the technique at the end of the proof of Lemma 14.21 It 
is a consequence of the variational property of X) namely that ||x||m ^ llV'llm with -0 a properly 
scaled bump function: t{j{a) := a | j . Computing its norm as we did in Lemma 14.21 gives 



<C2 



and the second inequality follows. 



□ 



We are now able to give the pointwise bound on the Lagrange function by applying Lemma [3. 81 (the 
zeros lemma). We note that when rM = oo, the first estimate in Proposition 14.51 holds throughout 
M (this is the case when M = M.'^). For the surface spline kernels the following proposition holds 
by applying [231 Corollary p. 130], while for Sobolev splines it follows in precisely the same way as 

for KM,m- 

Proposition 4.5. //H has mesh norm h < Fi/iq min(l, rM/2) then we have the estimate 



' h\ ^^'^/^ 
|x(a)| < Cm,M ( - 



exp 
exp 



/ dist(g,a) 

(?)] 



for dist(a, ^) < rj 
for dist(a, ^) > rj 



Proof. For a near to ^, say dist(^,a) < ytK^' apply Lemma \3M with r = "^YTh^ (i-e., on the 

ball b(a, 2p^), where the zeros of have mesh norm bounded by 2h) to achieve 



Fi/io 



m-d/2 



<Cr. 



1, 
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The second inequality is a consequence of Corollary 14.41 with T = 0. Thus, for a € b(,^, we 
have |x(a)| < C^.k 



m-d/2 



dist(^,a) > 



exp(-z/ — ^ 

A similar argument holds for < T = dist(a,^) < rM, where we can apply Lemma 13.81 (now 



with r 



to obtain 



lx(«)l < c, 



h 



Ti/io 



m~d/2 



IXlliy-(b(a,j^)) 



lxllw/-(b{a-,T-r^)C)- 



The second inequality follows because the ball b(a, is contained in b(,^, T — p^)'^- Applying 
Lemma 14.41 gives 



lx(a)| < C„ 



< a 



m-d/2 



exp 



h 



m-d/2 



exp 



/ dist(a, ^) 
V h 



For T > iM, we make use of the fact that b'^(a,rM) D h^{a,T), so we simply use the estimate for 
T = rM to obtain 

h \ ™~'^/^ 

exp 



\xia)\ < C 
since \\x\L,bHa,T) ^ \\x\L,bHa,VMy 



h 



□ 



4.3 The Lebesgue Constant is Bounded 



We are now in a position to prove our main theorem about the boundedness of the Lebesgue 
constant. It follows from the fast decay of the Lagrange functions by a standard argument that 
decomposes M en annuli, counting the elements of H in each annulus and balancing this against the 
influence of the far away Lagrange functions. Quasiuniformity is used in two ways. First, to obtain 
bounds on Lagrange functions that do not involve h/q, and, second, to estimate the cardinality of 
subsets of H. We note that f], a measurable subset of M, can be covered by a collection of small 
balls (b(^, g/2))^^„^^ that never overlap. Thus 

#(=nfi) < ,(b«,,/2)) < C-V"). (4.6) 

with C = Cm depending only on the manifold. 

Theorem 4.6 (Lebesgue Constant). Let M. be a complete, compact Riemannian manifold of di- 
mension d, and assume m > d/2. For a quasiuniform set H C M, with mesh ratio h/q < p, if 
\h < Ti/iQ min(l, rM/2), then the Lebesgue constant, L = sup^gj^ S^eH associated with 

i^m,Mi^^ bounded by a constant depending only on m, p and M. 
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Proof. Fix X. We first split the sum 



56S 
dist(^,o)<rM 



dist(5,a)>rM 



IX^HI =:/ + // 



into two parts. 



The "outer" sum, II, can be estimated by Holder's inequality using ()4.6p and Proposition 14.5 



// < Cq-'^n{M) exp 
which is bounded indepent of h. 



< Cp'^fi(M)h-^ exp 



The inner sum, I, is further decomposed into sums over sets Ek := H n a{a,k,h). By (j4.6p . 
the cardinality of Sfc is less than ^^(H n h{a,hk) < Cq~'^fi{h{a,hk)) < C{pk)'^. The Lagrange 
functions centered in are bounded by |x^(q^)I ^ C exp [— (A; — 1)] < Cexp[— z^A;], because 
dist(^, a) > {k — l)h. Thus / can be estimated by 



/ < E exp [-i^k] < Cp' 

k=i eeSfe 



oo 
k=l 



k'^ exp [-I'k] 



which is also bounded. 



□ 



Remark 4.7. We note that, when M = M'^, one may apply this theorem to the family of kernels 
known as the surface splines, as well as to the family known as the Sobolev (or Matern) splines. 
This follows in a straightforward way from Proposition 14. 5^ which can easily be shown to hold for 
the Lagrange functions associated to either of these kernels. The crucial difference is that, when 
M = W^, the radius of injectivity is rM = oo, and Lagrange function is bounded by ^ 



Cm,d \ 



m-d/2 



exp 



\(-a\ 



throughout 



It follows that sum in the preceding proof is 



controlled by the inner sum / only (and not //). 



5 Example: Interpolation on 

In this section, we provide the example of the kernel K2 §2 that gives rise to bounded Lebesgue 
constants for quasi-uniform data and whose approximation order can be given exactly. This 
will be accomplished by establishing that K2{x,y) =: (j){x ■ y) belongs to a class of SBFs studied 
in [24] whose L^, 1 < p < 00, approximation orders were explicitly obtained. For the readers 
convenience, we review relevant details from [24j which contains a full discussion of the material 
below. 

Let Hi denote the span of the spherical harmonics with fixed order i on S'^. The orthogonal 
projection Pi onto Hi is given by 

Pef = (/; Ye,m)ye,m- 
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Let Ld := a/A^ — Agd be the pseudo-differential operator where Agd is the Laplace-Beltrami oper- 



ator on S"' and = The Bessel potential Sobolev spaces have the norm given by 



oo 



£=0 LP(§d) 



Finally for /? = 1, 2, . . ., let Gp be the Green's function for L^, i.e., 

LdGnix-y) = 6y{x). 

More generally let (f)^ := Gp + Gp * tp where tp is an zonal function and let := span{(/)^(3; • 
^) : € S} where S has mesh ratio h/q < p. 

The following was given in \24:\ Theorem 6.8]. 

Theorem 5.1. Letl <p<oo, 7 > 0, f3 > 'j + d/p' 1/p' = 1). /// E H^, then with (j)^ 

and Ss as given above, and for quasiuniform centers H having mesh ratio p, 



We next wish to show that K2{x,y) = (j){x-y) defined on of the form cj) = ^4 + 04*7/), ip £ L^{S'^) 
and zonal. By the above theorem, <p will have L°° approximation order 4 for functions in / G C^(S^). 
For functions of lower smoothness, a simple IT- functional argument shows that for / G C^, 



To this end let (j) = (j){x, y) be the reproducing kernel for the Hilbert space with inner product 



dist^p(/,5H)<C/l''-Vll/ll 




{u,v)n := (V^n, V2i;)i2(52) + {Vu,Vv) ^2^32) + {u,v)l2^s^) 



where the Vu is the covariant derivative and 



u, the "Hessian," is given by 




Thus (n, (/))7v = u{x) where 



{u,(f>)N = (V^n, VV) + (Vn, V(^) + {u,(f>) 
= {u,L(t>) - {u,A(P) + {u,(P) 

= (-u, (L - A + =: {u,L(j)). 
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One next needs to identify L, which, in turn, wih lead to the identification of the Green's function 
(j). This requires calculating the Christoffel symbols. 

We will do this in (0, 92)-spherical coordinates, with 6 being the colatitude and (p is the azimuthal 
angle (longitude). The orthonormal basis for a tangent plane at (0,<^) comprises unit tangent 
vectors along and (p, {eg, e^}. With some work, it can be shown that 

V^u = uee eg®eg + esc 9{ug^ - cot 6){ee + (g) eg) 
( 1 

+ . 9 ^ u^^ + cot 

At this point we postulate that L is zonal in which case all partial derivatives of u with respect to 
ip are zero. It will be shown that L will have a zonal Green's function. In this case, V^u reduces to 



V^u = ugg eg (g) eg + cot Buge^ ® e^p 



V u ■ V V = uggvgg + CSC 9{sin 6 COS 9ug } [v^^ + COS 9 sin 9vg] 



and 
It follows that 

J V^U • V^W dfi = J UggVgg sm9 d9d(p 

+ J (csc^0 sin 9 cos 9ug){v^^ + sin 9 cos 9vg) sin 9 d9d'p. 
Upon integration by parts first in the (f variable and then in the 9 variable one obtains 



V sin 9d9d(p where 

Js2 y§2 



1 ( d^ f . d'^u\ d f , diji\ 
I^eu = — — - < -rr^ sm 9——^ — cot 9 cos 9— 



sm ( 



\ d02 \^ ^Q2 ) d9\ d9 J 



Thus the Green's function (p must be a solution to 

Lu={Lg-A + I)u (5.1) 
^ (sin 9u")" - (cos 9 cot 9u')' - (sin 9u')' + n = 



sm d sm t) sm ( 

for G (0,^). 

The equation (jS.ip has four analytic solutions in the interval (0,7r). At 9 = and 9 = it, it has 
regular singular points. We will use the method of Frobenius for analyzing the solution and the 
associated singularities [SJ p. 132]. 

In this method, a solution to (|5.ip is assumed to have a power series solution of the form 

oo 

u{9) = 6*^ ^ bk9'', where /? G C is unknown. 
fc=o 
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Substituting u{9) into (j5.ip results in the indicial equation, which is — 2)^ = in this case. 
Two solutions to ()5.ip are determined by the two roots /3 = and (3 = 2 and the corresponding 
equations for the hks. The other two of the four linearly independent solutions are then derived 
from these. Carrying out the details yields the four solutions below. 

ui{e) = pi{9^) and U2i9) = ln{0^)p2{0^), /3 = 0, (5.2) 
n^ie) = eWO^) and u^O) = 9^ (p^O^) + lnie^)p5{e^)) ,13 = 2. (5.3) 

The pj's are all power series convergent in a neighborhood of ^ = 0. In addition, they all satisfy 

pm + 0. 

It is important to note that these solutions have only even powers in their series. This is a 
consequence of the symmetry Q — > —Q of the differential equation and the method of Frobenius 
itself. 



There is another symmetry that preserves (j5.1|) : namely, Q — > t: — 6. This transforms the regular 
singular point 9 = n into the one at 9 = 0. Consequently, the four functions Uj{TT — 0), j = 1 to 4, 
are linearly independent solutions to (|5.ip valid near 9 = n. 

_ 4 
Thus the Green's function (p for L has the form (f)(9) = ajUj{9). Note that (p must have finite 

i=i 

norm, i.e., 

ml := (V2<^, V^,/.) + (V(A, V</)) + {(PA)- 

By Kondrakov's embedding theorem (see [2j), (p necessarily must be continuous (as a function of 
9) on [0, vr]. Thus (p cannot have either a U2{9) component or a U2{tt — 9) component. The singular 
behavior of the Green's function (p at 9 = is thus determined by U4{9). 

If there is any singular behavior at = vr, it will come from n4(7r — 9). This behavior is precisely the 
same type as at = 0. If the Green's function cp at n had the same singular behavior as at = 0, 
then L<p would exhibit nonzero distributional behavior at 9 = n, producing a second 5 function 
there. Since this doesn't occur for the reproducing kernel, the Green's function (p{9) is analytic 
near 9 = n. Here is a summary of what's been found: 

^(0) = / ^^'i^^') + ^^'^'3(^') + C9^ {{P4{e^) + H9^)P5i0^)), 9^0, 
[ Dp,{{7r - 9)^) + Ein - 9)^ps{{7r - 9)^), 9^n. 

We also mention that (p{9) is real analytic in 9 on the interval (0, tt]. (Of course, it is singular at 
9 = 0.) 

There is one final step that we need to take to determine the order of approximation; namely, we 
need to obtain the Green's function (p{9) in terms of t = cos9. Near 9 = 0, we write sm{9/2) = 
Y^(l — t)/2 and obtain 9 = 2 arcsin — t) /2, for which 9^{t) is analytic in a neighborhood of 
t = 1. At TT, we use tt — 9 = 2 arcsin {1 + 1)/2. This gives us (tt — 9{t))'^ analytic near t = —1. 
Inspecting the structure of the solution in ()5.4p shows that (p can be expressed in the following way 
as a function of t = cos 9: 

^{9{t)) = ao(l - t) ln(l - t) + • • • + aK(l - t)"" ln(l -t) + Mt), 
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with ao 7^ 0, and where, by adjusting K, fx can be taken to be as smooth as one wishes on [—1, 1]. 
To finish up, observe that in [24'j p. 8], this singularity corresponds to/3 = 2- l + 2 = 4. Thus the 
approximation rate is h'^. 
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